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Abstract. We consider the energy spectrum of a quasi-geostrophic model of 
forced, rotating turbulent flow. We provide a rigorous a priori bound E(k) < 
Ck~ 2 valid for wave numbers that are smaller than a wave number associated 
to the forcing injection scale. This upper bound separates this spectrum from 
the Kolmogorov-Kraichnan k~3 energy spectrum that is expected in a two- 
dimensional Navier-Stokes inverse cascade. Our bound provides theoretical 
support for the k~ 2 spectrum observed in recent experiments. 



The typical time scales associated with atmospheric flow over long dis- 
tances are much bigger than the time scales associated with Earth's rotation. 
This low Rossby number situation is characterized by a relative suppression 
of momentum transfer across vertical scales, and the organization of the 
flow in quasi-two dimensional strata. Cyclonic and anti-cyclonic vortical 
motion ensues in these layers, with dynamics in which strong, interacting 
vortices of many sizes are born, grow, and dissipate, over time scales that 
are long compared to the rotation time scale. The precise mathematical way 
of describing such a quasi-two dimensional picture is yet unclear. Energy 
spectra are some of the most robust quantitative indicators that one can 
use in order to distinguish between different classes of models. If a strictly 
two dimensional Navier-Stokes framework is adopted for rotating turbulence 
then the predicted energy spectra are a k~ 3 direct enstrophy cascade (at 
wave numbers larger than the wave number of the forces stirring the fluid) 
and a Kolmogorov-Kraichnan k~s inverse energy cascade spectrum at wave 
numbers that are smaller than the forcing wave numbers JIJ. Recent ex- 
periments |2j of rotating fluids find a different inverse energy cascade power 
spectrum: E ~ k~ 2 . This spectrum implies a steeper inverse energy cascade 
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than the one predicted by a strictly two-dimensional Kolmogorov-Kraichnan 
spectrum. The k~ 2 spectrum was observed for wave numbers of 10 _1 — 10° 
cm -1 . The experimental data showed the steeper k~ 2 spectrum clearly sep- 
arated from a k~3 spectrum that was fit to agree with E(ko) at the largest 
scale ko ~ lCT 1 cm -1 . 

The purpose of this letter is to describe a rigorous upper bound E(k) < 
Ck~ 2 valid in the inverse cascade region, in a quasi-geostrophic regime. 

The most important feature of strongly rotating fluids is the geostrophic 
balance between the Coriolis force and pressure gradients. This balance, 
valid only in a first approximation, imposes a two-dimensional time indepen- 
dent solution. The departure from this balance, to lowest order, has non- 
trivial dynamics and is described by quasi-geostrophic equations ||, which 
are quasi-two dimensional equations asserting the conservation of potential 
vorticity q subject to dynamical boundary conditions. The simplest of these 
can be written as 

d t 9 + v V6 + w E A6 = f (1) 

The two-dimensional velocity v is incompressible, V ■ v = 0. The dissipative 
term weA-9 has a coefficient we > that comes from Ekman pumping at 
the boundary. This coefficient has units of velocity and, for the situations 
we consider, where a non-trivial vertical flux is imposed, this coefficient is 
not vanishingly small. (In the experiment (0) the vertical velocities at the 
boundary are close in magnitude to the maximal measured velocities, and 
are in ranges of about 20-30 cm/s; the value for we is expected to match 
the same order of magnitude). The operator A can be described in Fourier 
representation as multiplication by the magnitude k of the wave number 
k — (^kx kyj j thcitj is! 

A0(k,t) =fc0(k,t). (2) 
The velocity is related to 9 by 

v(k,t) = (-k y ,k x )^9(k } t) (3) 
Thus, for each wave number 

|0(k,t)| = |v(k,t)|. (4) 
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The symbol / refers to the forcing term. This active scalar surface quasi- 
geostrophic equation has been studied both analytically and numerically || . 

We will analyze the energy spectrum using the Littlewood-Paley decom- 
position 0. The Littlewood-Paley decomposition is not orthogonal, but it 
is nearly so. Its use affords great flexibility in dealing with functions that 
involve many active scales: wave numbers are grouped in dyadic blocks and 
averages over the dyadic blocks are performed. The Littlewood-Paley de- 
composition is defined in terms of a smooth partition of unity in Fourier 
space. This partition is constructed starting from a nonnegative, nonin- 
creasing, radially symmetric function 0(o)(k) = 0( O )(fc), that equals 1 for 
k < !&o and vanishes for k > |& . The positive number k is just a 
reference wave number that fixes units. The argument of this template 
function is then dilated, setting <f>( m )(k) = 4>(2~ m k) and then the tem- 
plate is differenced, setting ipr \(k) = 4>(i){k) — <p^(k), and then ip( m )(k) = 
-0(o)(2~ m fc) = 4>( m +i){k) — 4>{ m )(k) for all integers m. The functions i/j^(k) 
are identically one for k G [|2 m fc , |2 m fc ], and vanish outside the inter- 
val [|2 m fc ,f2 m fc ]. The relationship 1 = (m) (k) + £^ =m Vn(k) holds for 
any integer m. One defines the Littlewood-Paley operators and A n 
as multiplication, in Fourier representation, by 0( m )(k) and, respectively by 
^(n)(k). Symbolically this means that the identity operator is written as 
I = S^ m > + Y^=m ^n- The Littlewood-Paley decomposition of a function 
F is F = S^F + Yl^m^-nF. For mean-zero functions that decay at in- 
finity, the terms S^ m ' becomes negligible when m — > — oo and therefore, for 
such functions one can write F = Yl™=-oc A n F. It is easy to see that for 
each fixed k > at most three A n do not vanish in their Fourier representa- 
tion at k (i.e. the conditions k G [|2 n /c 05 §2 n A<)] can be satisfied by at most 
three integers n, because n G [—1 + log 2 (^), 1 + log 2 (^-)]). The operators 
and A„ can be viewed as convolution operators. In particular, for ev- 
ery n = ±1, ±2, . . . , A n = f dh\l/( n )(h)5h- Here \l/( n ) is the function whose 
Fourier transform is ip( n ), ^( n ) = V'(n)' an d #h is the finite difference operator, 
(8\iF)(y) = F(r — h) — F(r). Thus A n is a weighted sum of finite difference 
operators at scale 2~ n k^ 1 in physical space, {k^ 1 provides thus an (arbitrary) 
length unit). 

The Littlewood-Paley decomposition of the solutions of the quasi-geostrophic 
equation is performed at each instance of time, 

oo 

6(T,t)= e in){*,t). 
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We wrote for ease of notation 6( n ) instead of A n #, thus, 



( n)M) = J dh¥ (n) (h)<f h (0)(r,t). 

The finite difference is taken at equal times, 6^(8) (r, t) = 8(r — h, t) — 6(r, t). 
There is an analogous decomposition for the velocity v and the forcing term 
/. In particular, in Fourier variables, the Littlewood-Paley components of v 
are given by vv^(k, t) = ■0( n )(k)v(k, t). The Littlewood-Paley spectrum || 
is 

E LP (k) = l Yl <I^(M)| 2 > (5) 



-l+log 2 (i)<n<l+log 2 (A) 
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where (...) is space-time average. The relation to the usual energy spectrum 
is straighforward. Because 

v(k,t)= J2 v w (k,t), (6) 

-l+log 2 (^)<n<l+log 2 (^) 

it follows that the usual energy spectrum 

E(k) = l(\v(k,t)\ 2 ) (7) 

satisfies 

E(k) < 3E LP (k). (8) 

Clearly, because the functions i])r n \ are non- negative and bounded by 1, one 
also has 

E(k) > l -E LP {k). (9) 
The temporal evolution of the system induces an evolution 

d t {n ) + v ■ V9 (n) + w E A6 {n) = R {n) (10) 

where 

R(n) = /(„) + J dh^ (n) (h)V h ■ (5 h (v)8 h 6) (11) 
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and fi n ) is the Littlewood-Paley component of the forcing term. Multiplying 
([T0| ) with 9( n ) and taking space-time average, one obtains the balance 



w E k( 9 {n) (k,t) > = (R (n) (T,t)6 (n) (T,t)). 
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Let us consider the case when the forcing term has a limited support in 
Fourier space, f(k, t) = for k ^ [k a ,kb], with < k a < kb < oo. The 
inverse cascade region will be described by wave numbers smaller than the 
minimal injection wave number k a . The inverse cascade region corresponds 
thus, in the Littlewood-Paley decomposition, to indices n > — oo that satisfy 
2 n+1 k < k a . We show now that the right hand side of the equation (|12|) is 
bounded above, uniformly for all such n > — oo. Because we are in a region 
where = 0, the term in the right-hand side of fll~2| ) can be written, after 
one integration by parts, as 



(R [n) {r,t)d {n) {r,t)) 



rfhV h ^ (n) (h) (5h(v) (r, t)5 h {6) (r, t)6 [n) (r, t)) . (13) 

This is a weighted sum of triple correlations. We will analyze each of the 
three terms involved in it, in an elementary but rigorous fashion. Because 
we aim at an upper bound, we will not try to optimize the prefactors. Using 
the Fourier inversion formula 

e (n) (r,t) = (2vr)- 2 J rfke jr -V (ri) (k)^(k,t), 

the term 9^ n ) is bounded pointwise by applying the Schwartz inequality: 

|^ (n) (r,t)|<(27r)- 2 ||^ (n) ||||^||, (14) 

with || • ■ ■ || the mean square norm. Using the fact that ipi n ) is a dilate of 

V>(o), we get 

|0(n)M)| <ci,2 n E{t)*, (15) 

where c\ = (2tt)- 2 / dk|^ (0) (k)| 2 and E(t) = j dr\6{r,t)\ 2 = / dr|v(r, t)| 2 
is the instaneous total energy. In the equality above we used Plancherel's 
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identity \\F\\ = (2tt) 1 \\F\\, and (Q). In order to bound the other two terms 
we note that, from Plancherel, we have 

J dr\5 h 6(r,t)\ 2 = (2n)- 2 J dk | e - lh ' k - 1| 2 | 0(k, t) \ 2 . 

Using |e"* h ' k — 1| 2 < 4hk, we deduce 

J dr\6 h 0(r,t)\ 2 <4h V (t) } (16) 

where r)(t) = J drO(r,t)A6(r,t). The term involving S^v is bounded using 
the same argument. In view of the bound is by the same quantity: 

J dr\S h v(r,t)\ 2 <Ah V (t), (17) 

Putting the three inequalities ([15], ^6), [l^) together with the Schwartz in- 
equality, we deduce that the triple correlation term that is integrated in ( 13 ) 
obeys 

\(S h (v)(r,t)8 h (9)(r,t)9 (n) (r,t))\ < 4c^2 n h(E(t^ V (t)) (18) 

In view of the fact that the functions \&( n ) are dilates of a fixed function, we 
deduce that 

\(R (n) (r,t)e (n) (r,t))\ < 2 n+1 C^. (19) 

Here E is the maximum total (not per unit volume) kinetic energy on the 
time interval, E = swp t E(t). The constant 

CV = 2c^ J dh/i|V h ^(o)| = c k (20) 

is proportional to ko and depends on the choice of the Littlewood-Paley 
template ?/>( ) only through the non-dimensional positive absolute constant 
Co- The number 

v = (v(t)) 

is related to the long time dissipation. It can be bound in terms of the forcing 
term using the total balance 

X d f 

--E(t) + w E r ] (t) = J dTf(T,t)6(r,t), 
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which follows from ([I]) after multiplication by 9 and integration. Writing the 
integral in Fourier variables and using the fact that the support of the forcing 
is bounded below by k a > one obtains the bound 

v^wfk-wf^m. (21) 

This bound diverges for very large scale forcing, i.e. when k a — > 0. Neverthe- 
less, because of the presence of the coefficient 2 ra+1 in ([19]) and the fact that 
2 n+1 ko < k a in the inverse cascade region, the total bound on the spectrum 
does not diverge as k a — ► 0: inserting ( |21| ) in ( |T9"D and using ( p0|) we get 



\(R (n) (T,t)6 {n) (r,t))\ < c w^Emf(r,t)\ 2 }. (22) 

Now, using ( p2|) in (|i~2"D and recalling the definition @ and the inequality (H) 
we obtain 

E(k) < Ok' 2 (23) 

for all k < k a . This is the main result of this letter. The constant has units 
of length per time squared and is given by 

C = 3c £^- 3 (|/(r,t)| 2 >. (24) 

The upper bound proved in this letter holds in greater generality than 
presented here. First of all, the spectrum of the forces need not be con- 
fined to the band [k a ,kb\. The role played by (ka)" 1 is then played by 
the ratio (k~ l \ f(k, t)| 2 ){(|/(k, t)| 2 )} -1 . Secondly, the results and methods 
apply to a much wider class of quasi-geostrophic equations. In fact, the 
quasi-geostrophic equation chosen here is the simplest version of a class of 
equations in which the potential vorticity q is advected by a a three dimen- 
sional velocity field that can be derived from a stream function ip (not to 
be confused with our Littlewood-Paley cutoff functions). The velocity has 
no vertical component v = (u,v) = (—d y ip,d x ip)- The potential vortic- 
ity q is advected following horizontal trajectories, d t q + v • Vg + (3v — F, 
where F includes sources and damping. The potential vorticity and the 
stream function are functions of three space variables (x,y,z); the relation 
Q = (d xx + d yy + d zz )ip closes this equation. The boundary conditions at z = 
are (P with 9 = d z ip. One decomposes ip — ips + i n a sum of a harmonic 
function ip B , (d xx + d yy + d zz )ipB = 0, <9 2 ?/>_b| 2= o = 6 and a function that sat- 
isfies homogeneous Neumann boundary conditions, (d xx + d yy + d zz )ipN = <?, 
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d z ipN\ z =o = 0- The ensuing equations on the boundary can be analyzed as 
above, using properties of the smooth evolution of q. 

Two main ingredients were used in the proof. The first one is the way 
in which the spectrum of the potential temperature is related to the energy 
spectrum. The second, and the essential ingredient is the fact that the re- 
laxation time at wave number k is roughly (w^fc) -1 in the range of wave 
numbers considered. This dependence is an important by-product of the 
quasi-geostrophic model. In contradistinction with direct cascade models 
where there is a dissipation anomaly, in the quasi-geostrophic models the 
coefficient we is not vanishingly small. Moreover, the physical forcing am- 
plitude obtained from an Ekman boundary layer is proportional to %. The 
explicit presence of the large scale term in the pref actor C is a reflection 
of the fact that the k~ 2 spectrum is modified near k = 0. The same is true for 
the Kolmogorov-Kraichnan spectrum: the integrals diverge and require an 
infrared cutoff. Such a cutoff can be achieved mathematically in two different 
ways. One may impose a smallest wave number k m i n and boundary condi- 
tions; or one can modifiy the dissipation law so that, in the limit k — > one 
has a finite relaxation time. In either case one can prove an a priori upper 
bound on E depending on the forcing and dissipation mechanism. The fact 
that the large scales are nearly conservative, with finite energy, was used in 
|2j to compare the k~ 2 and k~s spectra with the same largest scale energy. 
Our upper bound (|23|) confirms theoretically the separation of the two spec- 
tra when the injection length scales is small enough. This is indeed the case 
in the experiment J|: forcing was applied through 120 holes of diameters of 
.25 cm, some 8 times larger than the Ekman boundary layer length. 

In summary, we have proved that the energy spectrum of a forced surface 
quasi-geostrophic equation is bounded above by Ck~ 2 for wave numbers that 
are smaller than the force's injection wave number. Such a bound distin- 
guishes the quasi-geostrophic model from a two-dimensional Navier-Stokes 
model, and agrees with the recent experimental evidence of 0. 

The author acknowledges useful discussions with R. T. Pierrehumbert 
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